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ABSTRACT 

We present a comprehensive study of the properties of nearly coherent brightness oscillations (NCOs) 
in a large sample of type-I X-ray bursts observed from low mass X-ray binary 4U 1728-34. We have 
analyzed ~545 ks of data of this source obtained with the Rossi X-ray Timing Explorer over a 3 year 
span from 1996-1999. The data contain 38 bursts, 16 of which show NCOs. We find no NCOs present 
when the inferred mass accretion rate of the system is lowest. Furthermore, we define a measure of the 
strength of the oscillations and find that this strength increases with increasing mass accretion rate. This 
correlation is particularly evident within bursts detected only a few weeks apart, and becomes less clear 
for bursts separated in time by several months to years. The correlation we find for 4U 1728-34 between 
the NCOs and the inferred accretion rate of the system is similar to that found for KS 1731-260 by Muno 
et al., where the NCOs are only present at relatively high mass accretion rates. However, unlike the 
case for KS 1731-260, we find an anti-correlation in the bursts of 4U 1728-34 between the existence of 
episodes of photospheric radius expansion and the inferred mass accretion rate of the system. Moreover, 
we distinguish between NCOs present in the rise and in the decay phase of the burst and find that the 
bursts with NCOs only in the decay phase occur at intermediate accretion rates, while those with NCOs 
in both the rise and the decay phases are concentrated at high inferred mass accretion rates. We discuss 
these results in the context of the theory of thermonuclear bursts and propose intrinsic differences in the 
neutron stars in KS 1731-260 and 4U 1728-34 as the origin of the different behaviors. 

Subject headings: X-rays: bursts — stars: neutron — stars: individual (4U 1728-34) 



1. INTRODUCTION 

Low mass X-ray binaries (LMXBs), systems in which a 
compact object accretes matter via Roche lobe overfiow 
from a low mass (< 1 Mq) main sequence companion, 
provide a fascinating window into the physics and dy- 
namics of accretion disks and compact objects. When 
the compact object is a neutron star, these systems also 
permit detailed studies of the burning of the accreted ma- 
terial on the neutron star surface. Thermonuclear flashes 
(called type-I X-ray bursts) occur in many LMXBs and 
are caused by explosive unstable burning of the accreted 
matter. The properties of these type-I X-ray bursts are 
affected by the rate at which the matter is deposited onto 
the neutron star (see, e. g., Lewin, van Paradijs & Taam 
1993, for an extended review of X-ray bursts). Thus, 
in order to understand these bursts, it is necessary to 
investigate the effect that variations in the mass accre- 
tion rate (M) onto the neutron star have on them. The 
first attempts at characterizing M in LMXBs centered 
around the persistent flux of the sources, since intuitively 
the larger the amount of matter accreted, the greater the 
amount of X-ray flux emitted. However, it has become 
clear that the persistent flux in general may not always be 
a good measure of M (e. g., van der Klis 1995) and hence 
cannot be directly correlated with the burst properties. 

X-ray color-color diagrams (CDs) were introduced to 
study subtle spectral changes in LMXBs (see, e. g., van 

^ Also Goddard Space Flight Center, Code 662. 



der Klis 1995, and references therein). Hasinger & van 
der Klis (1989) studied the correlated rapid X-ray vari- 
ability and spectral behavior of LMXBs and showed that 
they naturally fall into two distinct categories. These two 
groups are called the atoll and Z sources, after the tracks 
they trace out in the CD. The motion of the sources along 
these tracks is thought to be caused by changes in the 
mass accretion rate (Hasinger & van der Klis 1989) and 
it has become clear that the position of the sources on 
these tracks is a better indicator of the M of the system 
than the persistent flux (see van der Klis 1995 for LMXBs 
in general but Mendez 2000 for 4U 1728-34 in particular). 
With the launch of the Rossi X-Ray Timing Explorer 
(RXTE) satellite, previously unknown high frequency 
timing phenomena were discovered and now offer us new 
views of LMXBs (see, e. g., van der Klis 2000, for an ex- 
tended review of the recently discovered timing phenom- 
ena in LMXBs). In particular, nearly coherent bright- 
ness oscillations (NCOs) with frequencies between 330 
and 589 Hz were discovered in Type I bursts from sev- 
eral sources (Strohmayer et al. 1996; Strohmayer et al. 
1997a; Smith, Morgan, & Bradt 1997; Zhang et al. 1997, 
1998; Markwardt, Strohmayer, & Swank 1999; Wijnands, 
Strohmayer & Franco 2000). These NCOs have extremely 
high coherence values [y j IS.v > 4000; Strohmayer & 
Markwardt 1999; Muno et al. 2000), high amplitudes 
(up to ~50% rms; Strohmayer, Zhang & Swank 1997b; 
Strohmayer et al. 1998b), and they exhibit long term 



frequency stability of up to years (Strohmayer, Zhang & 
Swank 1998a). These properties strongly argue for the 
rotation of the neutron star as the origin of the NCOs. 
In this interpretation, during a thermonuclear burst in- 
homogeneities in the fuel or in the fuel burning may arise 
causing "hot spots" on the surface of the neutron star. 
These hot spots rotate with the spin frequency of the 
star and periodically modulate the X-ray flux. 

These NCOs can be present in the rise of the bursts, 
very near the burst onset, and/or in the decay phase (or 
tail) of the bursts. When the NCOs are observed in the 
tail, their frequency can increase by 1 - 3 Hz, reach- 
ing an asymptotic limit as the bursts progress (see e. g., 
Strohmayer et al. 1998b). It has been proposed that this 
frequency increase is caused by thermal expansion of the 
thermonuclear burning layer on the stellar surface. Heat- 
ing from the burning expands the shell by several tens of 
meters at the onset of the burst. Due to angular momen- 
tum conservation this layer will slow down causing the 
observed frequency to decrease. During the decay of the 
burst, the burning layer settles back on the neutron star 
surface and spins up causing the frequencies to increase 
again. In this picture, the asymptotic limit represents the 
true neutron star spin frequency (Strohmayer, Zhang & 
Swank 1997b). The atoll source 4U 1636-53 has shown 
evidence for a subharmonic to the NCO frequency, indi- 
cating the possible existence of two, nearly antipodal hot 
spots (perhaps at the magnetic poles; Miller 1999). This 
detection of the subharmonic raises issues about how the 
burning front propagates on the surface of the neutron 
star and how the fuel can be concentrated in two antipo- 
dal regions. 

The burst oscillations in the atoll source KS 1731-260 
have been recently studied by Muno et al. (2000). They 
reported the presence of a correlation between the in- 
ferred M of the system and the presence of NCOs in the 
bursts. They find that the bursts from KS 1731-260 can 
be separated into two classes: the so called "fast" bursts, 
which occur at relatively high M, exhibit NCOs, and show 
episodes of photospheric radius expansion (PRE, the ex- 
pansion of the neutron star photosphere occurring when 
the burst luminosity is higher than the Eddington limit), 
and the "slow" bursts, which occur at lower M, do not 
exhibit NCOs, and do not show PRE episodes. Studies 
like this require a large sample of bursts observed over 
a considerable range of accretion rates. Such data bases 
have so far only been obtained for the persistent sources 
because the outbursts of the X-ray transients are rare 
and when they do occurr only a limited number of X-ray 
bursts are seen. In this paper, we will concentrate on the 
persistent source 4U 1728-34. 

The LMXB 4U 1728-34 is a prolific source of type I 
X-ray bursts (Basinska et al. 1984) and one of the best 
studied LMXBs. It has been classified as an atoll source 
(Hasinger & van der Klis 1989; Mendez & van der Klis 
1999) and it exhibits most of the behaviors seen in those 
systems. It was the first neutron star LMXB to show 
twin kilohertz quasi-periodic oscillations (kHz QPOs) 
in the persistent emission and the first to show NCOs 



(Strohmayer et al. 1996). The kHz QPOs in 4U 1728-34 
have been studied by several groups (Strohmayer et al. 
1996; Ford & van der KHs 1998; Mendez & van der 
Khs 1999; Di Salvo et al. 2000b) and the NCOs by 
Strohmayer et al. (1996), Strohmayer, Zhang & Swank 
(1997b, 1998a), and Strohmayer & Markwardt (1999). In 
this paper we present a study of the behavior of NCOs in 
a set of 38 type-I X-ray bursts from 4U 1728-34, observed 
over a wide range of mass accretion rates. We found that 
4U 1728-34 is similar to KS 1731-260 in that the NCOs 
are only present in bursts which occur at high M (Muno 
et al. 2000). However, contrary to the case of KS 1731- 
260, we find an anti-correlation between the presence of 
PRE episodes in the bursts and the inferred mass accre- 
tion rate of the system. We also find a more complicated 
behaviour with M for the NCOs in the rise and in the 
decay phases of bursts. 

This paper is organized as follows: in §2 we describe 
in detail the observations being used for this study. In 
§3 we explain the methods used to characterize the per- 
sistent emission and the burst properties, and the tech- 
niques used to study the oscillations. In §4 we discuss 
the spectral states (i. e., the CD) of the source. In §5 we 
discuss the behavior of the NCOs and study them in cor- 
relation with other burst characteristics and the spectral 
state of the source. Finally in §6 we discuss the implica- 
tions of these correlations for our current understanding 
of 4U 1728-34 in particular and for thermonuclear burn- 
ing on the surface of neutron stars in general. 

2. THE OBSERVATIONS 

The source 4U 1728-34 was observed with the Propor- 
tional Counter Array (PCA) on board RXTE on many 
occasions between 1996 February 15 and March 1 (here- 
after referred to as observation set AOl), 1997 Septem- 
ber 19 and October 1 (A02), and 1998 September 30 and 
1999 January 18 (A03). The actual observation dates are 
listed in Tables 1, 2 and 3, along with the main proper- 
ties of all 65 observations in our sample. In total ^^545 ks 
of on-source data were obtained. The AOl observations 
were obtained with epoch 1 gain settings for the PCA 
instrument, while A02 and A03 observations correspond 
to epoch 3 gain settings. Hereafter, when we refer to an 
observation we mean those data that can be identified 
by an unique observation identification number or ObsID 
(see Tables 1, 2, and 3). All 5 PCA detectors were on 
during all observations except for 3.6 ks during the ob- 
servation with ObsID 20083-01-01-020 and 3.6 ks during 
20083-01-03-020. No bursts were detected during those 
two intervals, so we have excluded those data from our 
study. The analysis of portions of our data have already 
been published in the literature. Strohmayer et al. (1996) 
reported the discovery of NCOs using burst 5 in our sam- 
ple and Strohmayer, Zhang & Swank (1997b) presented a 
temporal and spectral analysis of bursts 3, 4 and 5, and 
reported the existence of NCOs in bursts 4 and 5. Mendez 
& van der Klis (1999) published a CD of AOl and a por- 
tion of A02 data (see also Di Salvo et al. 2000b). 

A variety of data modes was available for each ob- 
servation set in addition to the Standard 1 and Stan- 



Table 1 
AOl Observations of 4U 1728-34 (RXTE PGA Epoch 1) 



Table 2 
A02 Observations of 4U 1728-34 (RXTE PGA Epoch 3) 





Date 


Start 


End 


rp a 

J- exp 




Burst 
















Obs. ID 
















(10073-01) 


1996 


Time 


Time 


(ks) 






Obs. ID 
(20083-01) 


Date 
1997 


Start 
Time 


End 
Time 


rp a 

-^exp 

(ks) 


^ avg 


Burst 


01-000 


Feb 15 


11:51 


18:49 


11.8 


4.1 


1 
















01-00 


Feb 15 


18:49 


00:02 


13.1 


4.2 


2 


01-00 


Sep 19 


05:52 


10:40 


11.3 


4.8 




02-000 


Feb 16 


00:02 


06:22 


13.1 


4.4 


3 


01-01 


Sep 19 


12:29 


15:02 


6.6 


4.8 


13 


02-00 


Feb 16 


06:22 


10:14 


10.1 


4.7 


4,5 


01-020 


Sep 20 


07:32 


14:03 


14.7 


4.1 


14 


03-000 


Feb 16 


15:48 


22:11 


12.3 


4.7 


6 


01-02 


Sep 20 


14:03 


16:29 


6.2 


4.8 




03-00 


Feb 16 


22:11 


01:19 


7.7 


4.8 




02-01 


Sep 21 


15:43 


21:05 


13.8 


5.0 


15,16 


04-000 


Feb 18 


11:08 


18:49 


13.4 


3.4 


7 


02-000 


Sep 22 


05:59 


13:59 


3.3 


4.4 


17 


04-00 


Feb 18 


18:49 


23:49 


10.5 


3.2 


8 


03-01 


Sep 23 


23:50 


00:41 


1.3 


4.8 




06-000 


Feb 22 


11:32 


18:53 


8.9 


2.0 




03-000 


Sep 24 


09:14 


17:14 


18.9 


5.7 




06-00 


Feb 22 


18:53 


00:02 


10.1 


3.3 


9 


03-00 


Sep 24 


17:14 


18:18 


3.7 


5.4 




07-000 


Feb 23 


21:16 


05:16 


18.2 


2.3 




04-00 


Sep 26 


12:29 


19:55 


17.9 


3.9 


18,19 


07-00 


Feb 24 


05:16 


06:45 


3.2 


2.3 


10 


04-01 


Sep 27 


09:18 


16:43 


18.1 


3.2 


20,21 


08-000 


Feb 24 


17:51 


01:26 


15.5 


2.4 


11 


03-020 


Sep 30 


04:33 


10:39 


13.2 


3.9 




08-00 


Feb 25 


01:26 


05:11 


8.5 


2.4 




03-02 


Sep 30 


10:39 


14:05 


7.1 


3.9 




09-000 


Feb 25 


20:45 


04:45 


17.9 


2.5 


12 


04-020 


Oct 01 


06:09 


14:09 


17.5 


3.6 




10-00 


Feb 29 


23:09 


00:02 


1.1 


2.8 




04-02 


Oct 01 


14:09 


15:18 


3.5 


3.7 




10-01 


Mar 01 


00:02 


05:45 


13.3 


2.8 




















" Total on-source time per i 


ObsID 









"■ Total on-source time per ObsID 

^ Average flux in units of 10~® erg s" 

tainties of 1% - 3% (90% confidence interval). 



^ Average fiux in units of 10 ^ erg s ^ cm ^ with uncer- 
tainties of 1% - 3% (90% confidence interval). 



dard 2 modes. Because the same data modes were 
not always available across observation sets, we de- 
scribe here the extant modes and how they were used 
to extract the quantities of interest to this study. The 
AOl data were obtained in one burst trigger mode 
(TLA_ls_10.249_ls_5000_F), two burst catcher modes 
(CB_125us_lM_0_249_H with 122 /is time resolution in 
one energy channel; CB_8ms_64M_0_249_II with 8 ms 
resolution in 64 channels) and one single bit mode 
(SB_125us_0_249_ls with 122 ^s resolution in one chan- 
nel). The burst catcher modes provide data typically 
for 1 s before the burst onset and for as long as the 
count rate remains above the trigger level. In our case 
the length of data ranges from 13 to 16 s for these 
modes. For A02 we had the same burst trigger and 
catcher modes available as during AOl, however, in- 
stead of the single bit mode we now had one event mode 
(E_16us_64M_0_ls with 16 /iS resolution and 64 channels) 
available. The data for A03 were obtained in three sin- 
gle bit modes (SB_125us_0_13_ls, SB_125us_14_17_ls and 
SB_125us_18_23_ls all with 122 ^s resolution and one 
channel) and one event mode (E_125us_64M_24_ls with 
122 /zs time resolution and 64 energy channels). When 
combined, all four data modes for A03 data cover the 
complete RXTE energy range. To generate the CD we 
extracted light curves from the Standard 2 data in 4 
energy ranges to calculate the colors. We also used the 



Standard 2 data to generate spectra of the persistent 
emission and extract average fluxes. The power den- 
sity spectra of the persistent emission were calculated 
from the single bit mode data for AOl, the event mode 
data for A02, and the data of the combined modes for 
A03. These same data modes were also used to search 
for NCOs in the bursts. Finally, the time series of the 
burst spectra were obtained from the event mode for 
A02, the combined data modes for A03, and the 8 ms 
resolution burst catcher mode for AOl (hence there is at 
most 16 s of time resolved spectra for each burst in AOl). 



3. THE ANALYSIS 

To characterize the persistent emission, we calculated 
average fluxes for each observation. The fluxes were ob- 
tained by fitting (using Xspec; Arnaud 1996) a cut-off 
power law spectrum to the background-subtracted Stan- 
dard 2 data. A Gaussian line around 6 keV was added 
in order to obtain acceptable fits. Since we are only in- 
terested in the fiux of the persistent emission, we did not 
investigate the significance of this feature further (but 
see, e. g., Di Salvo et al. 2000a). The resulting averaged 
fluxes in the 3-18 keV range are listed in Tables 1, 2, 
and 3. In order to verify the reliability of our fluxes we 
flt other spectral models (3-component models: a comp- 
tonized spectrum plus cither a multi-color disk blackbody 



Table 3 
A03 Observations of 4U 1728-34 (RXTE PCA Epoch 3) 



Obs. ID Date Start End Tea^p"" Favg^ Burst No. 
(30042-03) 1998-99 Time Time (ks) 



Table 4 
Properties of Bursts from 4U 1728-34 



Date 



Start Fpk'' -Fbb^ PRE rise tail 
Time UT NCO NCO 



01-00 


Sep 30 


07:11 


11:52 


10.4 


4.4 


22 


1 


1996 Feb 15 


17:58:11 


7.3 


8.9(8) 


Y 




Y 


02-00 


Oct 07 


05:05 


09:49 


11.1 


5.9 




2 


1996 Feb 15 


21:10:21 


7.9 


9.4(8) 


Y 


- 


Y 


01-01 


Oct 23 


22:35 


23:08 


1.9 


3.8 




3 


1996 Feb 16 


03:57:10 


4.5 


5.4(6) 


- 


- 


- 


02-01 


Oct 24 


00:15 


01:09 


2.5 


3.8 




4 


1996 Feb 16 


06:51:08 


7.2 


11.(1) 


- 


Y 


Y 


03-01 


Oct 24 


01:47 


03:46 


4.8 


3.8 


23 


5 


1996 Feb 16 


10:00:44 


7.2 


11.0(7) 


- 


Y 


Y 


01-02 


Oct 26 


20:57 


21:34 


1.7 


3.9 




6 


1996 Feb 16 


19:27:11 


7.8 


10.(1) 


- 


- 


Y 


01-03 


Oct 26 


22:33 


23:10 


2.1 


3.9 




7 


1996 Feb 18 


17:31:50 


7.7 


11.(1) 


Y 


- 


Y 


04-00 


Oct 27 


00:09 


03:19 


6.9 


3.9 




8 


1996 Feb 18 


21:28:50 


8.1 


12.(1) 


Y 


- 


- 


01-04 


Oct 29 


01:44 


02:43 


3.5 


4.0 




9 


1996 Feb 22 


23:09:12 


9.2 


14.(1) 


Y 


- 


- 


05-00 


Oct 29 


03:20 


06:37 


7.0 


3.9 




10 


1996 Feb 24 


05:46:24 


8.7 


14.(1) 


Y 


- 


- 


06-00 


Nov 01 


02:30 


07:26 


9.8 


4.0 


24,25 


11 


1996 Feb 24 


17:51:49 


8.8 


13.(1) 


Y 


- 


- 


07-01 


Nov 02 


22:29 


23:30 


3.6 


4.0 


26 


12 


1996 Feb 25 


23:17:07 


8.8 


13.(1) 


Y 


- 


- 


07-00 


Nov 03 


00:05 


03:18 


7.4 


4.0 


27 


















08-00 


Nov 05 


04:00 


07:29 


6.7 


4.1 




13 


1997 Sep 19 


12:32:56 


5.3 


6.4(7) 


- 


Y 


- 


09-00 


Nov 10 


13:17 


13:55 


2.2 


4.4 




14 


1997 Sep 20 


10:08:50 


3.4 


4.0(5) 


- 


Y 


Y 


09-02 


Nov 10 


16:31 


17:07 


2.1 


4.4 




15 


1997 Sep 21 


15:45:28 


4.8 


5.7(6) 


- 


Y 


- 


10-00 


Nov 10 


18:26 


20:19 


4.9 


4.4 


28 


16 


1997 Sep 21 


18:11:04 


4.6 


5.5(6) 


- 


Y 


Y 


10-01 


Nov 10 


21:07 


21:54 


2.8 


4.4 




17 


1997 Sep 22 


06:42:51 


3.0 


3.6(5) 


- 


Y 


Y 


11-00 


Nov 10 


22:26 


03:16 


10.8 


4.3 


29 


18 


1997 Sep 26 


14:44:09 


4.9 


5.8(7) 


- 


Y 


Y 


12-00 


Nov 11 


14:50 


21:52 


17.9 


4.3 


30,31 


19 


1997 Sep 26 


17:29:47 


4.7 


5.5(6) 


- 


Y 


Y 


13-00 


Nov 11 


22:26 


03:13 


10.4 


4.4 


32 


20 


1997 Sep 27 


11:19:06 


8.5 


10.2(8) 


Y 


- 


Y 


14-02 


Nov 16 


06:34 


07:00 


1.5 


5.2 




21 


1997 Sep 27 


15:54:03 


7.9 


9.4(8) 


Y 


- 


Y 


14-01 


Nov 16 


08:14 


09:08 


3.2 


5.2 




















14-00 


Nov 16 


09:52 


10:47 


3.2 


5.1 


33 


22 


1998 Sep 30 


10:08:54 


9.4 


11.(1) 


Y 


- 


- 


15-00 


Nov 16 


13:09 


19:15 


14.9 


5.1 


34 


23 


1998 Oct 24 


02:32:28 


9.3 


12.(1) 


- 


- 


- 


16-00 


Nov 17 


00:12 


07:36 


17.1 


4.5 


35 


24 


1998 Nov 01 


03:47:55 


9.8 


12.(1) 


- 


- 


- 


17-00 


Nov 17 


13:02 


14:00 


3.4 


4.0 


36 


25 


1998 Nov 01 


06:59:11 


9.9 


12.(1) 


Y 


- 


- 


18-00 


Jan 16 


07:51 


11:52 


9.9 


3.1 




26 


1998 Nov 02 


22:42:59 


9.7 


12.(1) 


Y 


- 


- 


19-01 


Jan 17 


04:50 


05:34 


2.6 


3.1 




27 


1998 Nov 03 


02:02:24 


9.5 


11.(1) 


Y 


- 


- 


19-00 


Jan 17 


06:11 


10:10 


10.0 


3.2 


37 


28 


1998 Nov 10 


19:33:39 


8.3 


10.(1) 


- 


- 


- 


20-00 


Jan 18 


22:39 


01:18 


5.4 


3.3 


38 


29 


1998 Nov 11 


00:37:33 


9.2 


11.(1) 


Y 


- 


- 
















30 
31 
32 
33 
34 


1998 Nov 11 
1998 Nov 11 
1998 Nov 12 
1998 Nov 16 
1998 Nov 16 


16:20:37 
19:32:04 
00:19:54 
10:10:25 
16:09:03 


8.9 
8.8 
8.6 
8.4 
8.6 


11.(1) 
11.(1) 
10.(1) 
10.(1) 
10.(1) 


Y 
Y 

Y 
Y 


- 


- 


"■ Total on-source time per ObsID 

^ Average flux in units of 10~^ erg s^^ cm^^ 

tainties of 1% - 3% (90% confidence interva 


with 


uncer- 


Y 
















35 


1998 Nov 17 


07:30:33 


8.0 


10.(1) 


- 


- 


- 
















36 


1998 Nov 17 


13:44:06 


8.8 


11.(1) 


Y 


- 


- 
















37 


1999 Jan 17 


06:12:01 


8.6 


10.(1) 


Y 


- 


- 
















38 


1999 Jan 18 


23:46:33 


9.6 


12.(1) 


Y 


- 


- 



° Measured peak flux in the 3—18 kcV range in units of 
10~^ erg s~ cm~^. 

'' Peak bolometric Hux from best-fit black body model, 
in units of 10^* erg s^^ cm^^. 



or a plain black body plus a Gaussian line around 6 keV; 
see, e. g., Di Salvo et al. 2000a). We find that these alter- 
native models give the same average flux of the persistent 
emission within the quoted errors. 



To generate the CD, for AOl we define the soft color 
as the logarithm of the ratio of the count rate in the en- 
ergy range 2.9-4.5 keV to the count rate in the range 
4.5-6.3 keV. Similarly we define the hard color as the 




perform a standard spectrum analysis of these bursts us- 
ing Xspec. We find that the spectra of X-ray bursts are 
well fitted by a simple black body model (Swank et al. 
1977). As a background spectrum we used the spectrum 
of an interval of 10 s of the persistent emission before 
the burst onset. ^ While the true background spectrum 
during a burst may differ from that of the persistent emis- 
sion before the burst (van Paradijs & Lewin 1986) this 
procedure will suffice to establish the existence of PRE 
episodes. 

From the spectral fits we extracted time series of tem- 
perature, radius of the emitting region, and measured 
flux. The peak measured fluxes are listed in Table 4. 
From the best fit parameters and the black body as- 
sumption we can calculate a peak bolometric flux using 
Fbb — A-kR^gT^, which we also list in Table 4. We veri- 
fied the existence of PRE episodes by calculating hardness 
ratios as a function of time during the burst. In Figure 1 
we show a typical burst with a PRE episode. The drop 
in the hardness ratio in (6) indicates the beginning of the 
radius expansion phase. This behaviour is confirmed by a 
simultaneous decrease in the black body temperature (c) 
and increase in the inferred radius of the emitting region 
(d) while the bolometric flux remains constant (e). 

To search for oscillations in the bursts, we employ 
the Z^ statistic described by Strohmayer & Markwardt 
(1999, see also Buccheri et al. 1983). The main reasons 
for this choice are that it does not require binning and 
its probability density is x^ distributed even for small 
numbers of events (X-ray photons in our case). Since 
the oscillations are highly sinusoidal, we consider only Z^ 
(see Strohmayer & Markwardt 1999). We generate data 
segments of 2 s in duration with centroids separated by 
0.25 s (hence they are not independent), and search for 
power in the vicinity of the previously found oscillation 
frequency of 363 Hz (Strohmayer et al. 1996). We set the 
detection threshold at Z^ = 14. This is equivalent to a 

single trial significance of a little better than 1 x 10^'^. 
Since pulsations are known to exist in this source we do 
not need to search a wide frequency space and we there- 
fore have few effective trials. Using this criteria we found 
that 16 out of our sample of 38 bursts exhibit NCOs. 
Typical examples of bursts which show NCOs are shown 
in Figure 2. We established a measure of the strength of 
the oscillations by integrating the RMS amplitude of the 
oscillations over their duration in the burst, ignoring the 
wings of the power features which fall below our detection 
threshold. 

4. THE COLOR-COLOR DIAGRAM OF 4U 1728-34 

The CD of 4U 1728-34 in Figure 3 a shows a track typ- 
ical of atoll sources (Hasinger & van der Klis 1989, see 
also Mendez 2000). It is important to note that the shift 
in the CD of the track for AOl with respect to the tracks 
of A02 and A03 is not an intrinsic behavior of the source 
but due to the difference in gain settings. Our observa- 

^ Except for AOl observations where the only data mode available with high enough time and energy resolution was a burst catcher mode 
with only up to 16 s per burst. 

^ Except again for AOl observations where only up to 2 s of data was available prior to burst onset. 
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Fig. 1. — Typical burst (burst 22) with photospheric radius ex- 
pansion, (a) The count rate (in units of 10'*), (b) the logarithm of 
the ratio of counts in the 9-18 keV range to that in the 3-4 keV 
range, (c) the black body temperature, (d) the inferred radius of 
the emitting region, and (e) the bolometric flux (in units of 10^** 
erg cm~^ s~^) as a function of the time in the burst. The dashed 
vertical line marks the time at which the burst reached peak count 
rate. 



logarithm of the ratio of the count rate in the range 6.3- 
9.0 keV to that in the range 9.0-18.2 keV. For A02 and 
A03 the energy ranges are slightly different: 3.0-4.4 keV 
to 4.4-6.2 keV for the soft color and 6.2-9.1 keV to 9.1- 
18.2 keV for the hard color. The colors were calculated 
from background-subtracted light curves from the Stan- 
dard 2 data and hence have 16 s time resolution. How- 
ever, to reduce the errors on the colors we averaged them 
in 256 s intervals. 

We searched for bursts by generating light curves of 
the Standard 1 data with 0.125 s resolution. We found 
38 bursts in our data set. We have numbered the bursts 
according to time and have listed their main properties in 
Table 4. The start time of the burst is determined from 
the first time bin in the background-subtracted light curve 
with count rate above 1000 counts s~^. We characterize 
the bursts by obtaining time resolved series of spectra. 
We extracted dead-time corrected spectra of 0.125 s du- 
ration for a time interval of 20 s during the bursts.^ We 
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Fig. 2. — Representative dynamical power spectra for bursts with oscillations only in the rise of the burst (left), in both the rise and the 
decay portions of the burst (middle) and only in the decay phase of the burst (right). The first contour plotted corresponds to Z^ = 14. 
The contours increase in steps of 2 until Z^ = 20 and then they increase in steps of 5 until Z^ = 65. In the figures, the burst profiles are 
overplotted in the same arbitrary units. 



tions of 4U 1728-34 sampled both the island state and the 
banana branch, though not during a single observation 
set. Both AOl and A03 exhibit the complete island state 
track, from the extreme island state at the upper right of 
the CD, believed to be the state with the lowest mass ac- 
cretion rate, to just before the island state-banana branch 
vertex (hereafter simply referred to as 'the vertex'). Only 
during our A02 observations was 4U 1728-34 found in 
the banana branch. This observation set covers the en- 
tire banana branch, from the lower banana branch on the 
left corner to the upper banana branch on the bottom 
rightmost portion of the track, where the mass accretion 
rate is thought to be the highest. Our CD looks very simi- 
lar in shape to those published previously for 4U 1728-34 
(Mendez & van der Khs 1999; Piraino et al. 2000; Di 
Salvo et al. 20G0b). Our island state appears somewhat 
longer, but this difference is due to the different choice of 
colors and disappears when similar colors are used. 

To determine the position of the source in the CD prior 
to each burst, we determine the color of the source from 
an interval of 320 s before the burst, ending 16 s prior to 
the onset of the burst. The results can be seen in Fig- 
ure 3 a where the numbers indicate the spectral state of 
the source just prior to the occurrence of each burst. From 
the figure it is clear that we detected bursts throughout 
almost the entire track, from the extreme island state to 
the upper banana branch although not in the uppermost 
part of the upper banana branch. However, it is not yet 
clear whether this is due to a real lack of bursts at the 
highest M or due to the limited amount of data so far 
obtained when 4U 1728-34 was at such high accretion 
rates. 

To verify the spectral state of the source prior to each 
burst, average power density spectra (PDS) were obtained 
by integrating over the entire length of the continuous 
data train prior to each burst. Figure 4a shows the PDS 



of the source prior to burst 22, which is typical 
of the source when it is in the extreme island state. 
This PDS clearly shows band-limited noise at frequencies 
above 100 Hz, with no high frequency QPOs, a bump 
around 10 Hz, a flat spectrum at low frequencies, and a 
break around 1 Hz, typical for the island state (see, e. g., 
van der Klis 1995). At the vertex, e. g., prior to burst 21, 
the source PDS (Fig. Ah) shows a kHz QPO and a break 
around 20 Hz typical of the lower banana branch in tran- 
sition to the island state. Figure 4c shows the PDS for 
the source prior to burst 15, when 4U 1728 34 is in the 
upper banana branch. The PDS shows a typical feature- 
less power law noise component. As mentioned above, the 
difference in gain settings between epoch 1 and 3 causes 
a shift in the island state tracks between AOl and A03. 
To see where the AOl track might join the A02 track if 
there were no gain differences, we show in Figure Ad the 
PDS of the source prior to burst 3. This PDS is very sim- 
ilar to that of burst 21 (Fig. 46), indicating a transition 
between the island state and the lower banana branch. 

5. THE BEHAVIOR OF BURST OSCILLATIONS 
5.1. The NCOs and the Inferred Mass Accretion Rate 

As shown in Figure 5, we find no correlation between 
the average flux of the persistent emission prior to the 
bursts and the existence or the strength of the NCOs. 
However, if the presence and possibly the strength of the 
NCOs is correlated with the mass accretion rate, as was 
the case in KS 1731-260 (Muno et al. 2000), then the 
lack of correlation is not surprising, because, as already 
mentioned in § 1, the flux is not a robust tracer of the 
mass accretion rate. Instead, a better measure of M is 
the position of the source on the atoll track in the CD. 

In Figure 36, we plot all bursts in the CD; with dia- 
monds the bursts that do not show NCOs, and with cir- 
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Fig. 3. — a) Color-color diagram for 4U 1728-34 showing a typical atoll track and counterclockwise from upper right, the extreme island 
state (EIS), the island state (IS), the lower banana branch (LB), and the upper banana branch (UB). The stars correspond to data obtained 
with epoch 1 instrument gain settings (AOl data), while the exes (A02) and crosses (A03) belong to epoch 2. The shift in the tracks 
between AOl and A03 is due to these gain differences and is not intrinsic to the source. Each point represents 256 s of data. The numbers 
indicate the state of the source prior to the corresponding burst, b) The same diagram now showing the presence and the strength of the 
nearly coherent brightness oscillations (NCOs) in each burst. The diamonds are bursts without NCOs. The circles correspond to bursts 
with NCOs and the size of the circle is proportional to the strength of the NCOs (see text), c) The same diagram showing the presence of 
episodes of photospheric radius expansion (PRE). Bursts that exhibit PRE episodes are shown with squares, while those that do not are 
shown with diamonds, d) The same diagram showing the presence of the NCOs only in the cooling tail of the bursts (squares), both in the 
tail and the rise of the bursts (triangles), and only in the rise of the bursts (circles). Diamonds mark the bursts without detectable NCOs. 



cles those that do. For bursts with osciUations the size of 
the cireles is proportional to the strength of their NCOs 
as defined in § 3. From this figure it is clear that bursts 
with NCOs occur only on the bottom part of the dia- 
gram where the accretion rate is thought to be highest, 
i. e., near the vertex and on the banana branch. Close 
inspection of this figure reveals that the strength of the 



NCOs generally increases with decreasing hard color, i. e., 
the NCOs become stronger as the inferred mass accretion 
rate of the system increases. This relation is especially ev- 
ident for bursts within a single observation set, and hence 
separated in time by only a few days to a few weeks, but 
it is not so clear for bursts from different observation sets, 
separated in time by up to years (i. e., between the AOl 
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Fig. 4. — Representative power density spectra from bursts occurring when 4U 1728-34 was in the extreme island state (a), vertex 
between island state and banana branch (b, d), and upper banana branch (c). 



and the A02 data sets). 

5.2. The NCOs and the Burst Flux 

In Figure 6 we show the peak burst fluxes for the bursts 
in our sample. Bursts without NCOs are marked with 
asterisks, while bursts with NCOs are shown with circles 
whose radius is proportional to the strength of the oscilla- 
tions. Burst 3 which shows an unusual count rate profile 
is marked with a square. The first trend evident from Fig- 
ure 6 is that the bursts in our A02 observations (except 
for bursts 20 and 21) have considerably lower peak fluxes 
than those in our AOl and A03 observations. We also 
observe that the bursts without NCOs have on average 
higher peak fluxes. Furthermore, among the bursts with 
oscillations in a single observation set (AOl or A02), the 
dimmest bursts also have the strongest oscillations. This 
correlation does not hold when we consider bursts from 
AOl and A02 together. 

5.3. The NCOs and Bursts with Photospheric Radius 
Expansion 

In their study of KS 1731-260, Muno et al. (2000) found 
that all the bursts with PRE episodes occur at high mass 
accretion rates and they all show NCOs. They go fur- 
ther to distinguish between oscillations occurring during 
the rise of the burst and those occurring during the cool- 
ing tail of the burst. They note that oscillations in the 



rise of bursts are not correlated with PRE episodes since 
burst 7 in their sample exhibits such NCOs but not a PRE 
episode. However, their data suggests that all bursts with 
NCOs in the decay portion of the burst do exhibit PRE 
episodes. 

A careful look at Figure 3 c reveals a different picture 
for 4U 1728-34. Bursts with PRE episodes (shown with 
squares) occur throughout the entire island state and at 
the vertex, but not in the banana branch itself. This 
is precisely the opposite of the behavior in KS 1731-260 
where all bursts with PRE episodes lie in the banana 
branch, at the highest M. We also note that all 7 bursts 
from our sample that lie in the banana branch (but not 
bursts 20 and 21 which lie near the vertex) have very low 
peak measured fluxes (see Fig. 6) while all the PRE bursts 
in the KS 1731-260 study are "fast" bursts which show 
high peak fluxes. 

The behavior of the bursts in 4U 1728-34 is richer still. 
Muno et al. (2000) hinted at a possible difference between 
the NCOs observed during the rise of bursts and those ob- 
served in the cooling tails of the bursts. If we make this 
distinction we find in our sample 2 bursts with NCOs 
only in the rise, 7 bursts with NCOs in both the rise and 
the tail, and 7 bursts with NCOs only in the tail. Rep- 
resentative bursts of each of these categories are shown 
in Figure 2. In Figure 3d we can see that these distinc- 
tions in the NCOs are also borne out on the CD. Bursts 
with NCOs only in the decay tail of the burst are seen 
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Fig. 5. — Average persistent emission fluxes obtained from a power 
law fit to the data prior to each burst. The dashed lines separate the 
bursts during A02 observations from those of AOl and A03. Bursts 
without NCOs are shown with asterisks, while bursts with NCOs are 
shown with eircles whose size is proportional to the strength of the 
oscillations (see text). Burst 3 which exhibits an unusual light curve 
is marked with a square. 
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Fig. 6. — Peak burst fluxes obtained from a black body fit to 
running spectra of 0.125 s duration. The dashed lines separate the 
bursts during A02 observations from those of AOl and A03. Bursts 
without NCOs are shown with asterisks, while bursts with NCOs are 
shown with circles whose size is proportional to the strength of the 
oscillations (see text). Burst 3 which exhibits an unusual light curve 
is marked with a square. 



(marked with squares) at intermediate M, in the island 
state and the vertex. Those bursts with NCOs in both 
the rise and the tail of the bursts are found (marked with 
triangles) only at high M, in the banana branch. The 
two bursts with NCOs only in the rise of the burst are 
shown with circles and they appear in the upper banana 
branch, at the highest inferred M . However, with only 
two such bursts it is difficult to determine whether or not 
their location only in the UB branch is significant. The 
bursts without NCOs are shown with diamonds. 

6. DISCUSSION 

We have searched three years of archival RXTE data of 
the neutron star LMXB and atoll source 4U 1728-34 for 
type-I X-ray bursts. We have found 38 bursts, 16 of which 
show nearly coherent brightness oscillations (NCOs) . The 
presence or absence of the NCOs in the bursts is not well 
correlated with the persistent flux of the source, but it 
is very well correlated with the position of the source on 
the atoll track in the CD. The NCOs are only present in 
the bursts that occur when the source is near the island 
state-banana branch vertex or on the banana branch. It is 
thought (Hasinger & van der Klis 1989) that the position 
of the source in the atoll track is determined by the mass 
accretion rate. If true, then the presence of the NCOs is 
very well correlated with the inferred mass accretion rate 
and the lack of any correlation with the persistent X-ray 
flux is further proof that the persistent flux is not a good 
measure of the accretion rate of the system. 

A simflar study for KS 1731-260 by Muno et al. (2000), 
although with less bursts (9 bursts, 5 with NCOs), also 
found that the NCOs are only present in bursts which 



occur at relatively high mass accretion rates. However, 
Muno et al. (2000) also report that the bursts of KS 1731- 
260 which have episodes of photospheric radius expan- 
sion (PRE) occur at these high accretion rates and that 
they all exhibit NCOs. Such a correlation can be ex- 
cluded for 4U 1728-34. When we considered bursts that 
exhibit (PRE) episodes, we discovered a strong anti- 
correlation with inferred mass accretion rate: bursts with 
PRE episodes occur only when the source is at its low- 
est and intermediate mass accretion rates and not at its 
highest accretion rates. Also, the presence of NCOs dur- 
ing the bursts is not correlated with the presence of PRE 
episodes: both bursts with and without NCOs can exhibit 
periods of PRE. Yet another difference we found is that 
the bursts which occured during relatively high mass ac- 
cretion rates were the dimmest in our sample, while the 
bursts in KS 1731-260 were bright when the accretion 
rate in that system was high (Muno et al. 2000). It is 
unclear why the positions of bursts in the CD are very 
similar with respect to the existence of NCOs in both 
4U 1728-34 and KS 1731-260, but other burst properties 
differ so signiflcantly with CD position (and thus with the 
inferred mass accretion rate). 

One clue to these differences may be found in the the- 
ory of thermonuclear bursts. Atoll sources accrete mass 
at rates in the range 0.01 MEdd< Af < 0.1 MEdd (Hanawa 
& Fujimoto 1982). In this regime, the model of Fujimoto, 
Hanawa & Miyaji (1981) for thermonuclear bursts (see 
also Bildsten 1997) predicts two main behaviors domi- 
nated by the accretion rate. At relatively low mass ac- 
cretion rates, the hydrogen on the neutron star surface 
burns stably forming a layer of pure helium. As more 
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matter accumulates this layer is compressed and heated 
until it ignites unstably, producing a pure helium burst. 
As the mass accretion rate increases, hydrogen is accreted 
faster than it can be burned and when the helium flash 
ignities it does so in a mixed hydrogen/helium environ- 
ment. For atoll sources, bursts occurring in the island 
state at low inferred mass accretion rates would then be 
the result of pure helium flashes, with the fraction of hy- 
drogen fuel increasing as the source progresses toward the 
banana branch. 

Because of the relative timescales for energy release of 
the reactions involved (the CNO cycle for hydrogen burn- 
ing is limited by the rate of /3-decay and thus much slower 
than the triple alpha process for helium burning), bursts 
with high fractions of helium are expected to release their 
energy faster and reach higher peak fluxes. This behav- 
ior seems borne out in the bursts from 4U 1728-34. From 
Table 4 and Figure 3 a we can see that in our AOl ob- 
servations, the bursts in the extreme island state (bursts 
9-12) have higher peak fluxes than the rest of the bursts 
in that observation set. Similarly within the A03 obser- 
vations the bursts at the lowest M (22-27 and 38) also 
have the highest peak fluxes among their group. It is nat- 
ural then that episodes of photospheric radius expansion 
are anti-correlated with mass accretion rate; the lower the 
M, the higher the helium fraction, the faster the energy 
release and the brighter the bursts, making it more likely 
that they will reach the Eddington limit and expand the 
photosphere. 

In the case of KS 1731-260, Muno et al. (2000) argue 
that the inverse behavior of this source with respect to 
the theory's predictions may be explained by consider- 
ing not the overall mass accretion rate but the rate of 
mass accretion local to the fuel surface. If the area over 
which the matter is being accreted increases with mass 
accretion rate, the local accretion rate may in fact de- 
crease, causing the inverse behavior (Bildsten 2000). The 
question then becomes, why does the area over which 
the accreted matter is distributed increase with increased 
M in KS 1731-260 but not in 4U 1728-34? The answer 
may lie in the intrinsic differences of the neutron stars 
in each of these systems. For example, a significant dif- 
ference in the magnetic field (i. e., its strength and ori- 
entation) of the neutron star in 4U 1728-34 compared 
to that in KS 1731-260 may inhibit the motion of mat- 
ter from the equator to higher latitudes. Alternatively, 
if the frequency of the NCOs is indeed the true spin fre- 
quency of the neutron star, then KS 1731-260 with an 
NCO asymptotic frequency of 524 Hz is spinning nearly 
45% faster than 4U 1728-34 (with an asymptotic NCO 
frequency of 363 Hz) . This may reduce the effective grav- 
ity at the equator of the neutron star in KS 1731-260 and 
allow more significant spread of the accreted matter. 

If KS 1731-260 and 4U 1728-34 are so different with 
respect to their burst behavior, why then are burst oscil- 
lations not observed at low inferred mass accretion rates 
in cither of these sources? Evidently the mechanism pro- 
ducing the oscillations is strongly influenced by the mass 
accretion rate of the system. Since the NCOs are likely 



closely tied to the properties of the nuclear burning it 
is very possible that we are seeing the direct influence 
of mass accretion rate on the burning behavior. Muno 
et al. (2000) point to recent work by Cumming & Bild- 
sten (2000) which indicates that modulations of the burst 
emission are more likely observable during helium-rich 
bursts, which they argue occur in the banana branch for 
KS 1731-260 where the NCOs are observed in this source. 
This cannot be the case for 4U 1728-34 as we have argued 
that the bursts with NCOs in 4U 1728-34 have in fact 
high fractions of hydrogen mixed in which causes them 
to be dim. Further work on possible "hot spot" or other 
alternative mechanisms will greatly illuminate this issue. 

Although it is likely that a low mass accretion rate 
directly inhibits the oscillation production mechanism, 
we explore the possibility that these oscillations are in 
fact present during the burst but that environmental 
changes (due to the differences in M) cause the oscil- 
lations to become undetectable at low M . Recently, the 
X-ray eclipsing source X 1658-298 was found to also ex- 
hibit NCOs during five bursts (Wijnands, Strohmayer & 
Franco 2000). Several more bursts were found which did 
not show these oscillations. However, there is strong ev- 
idence that the absence of NCOs in those bursts was 
caused by environmental changes in the system. In 
X 1658-298, the oscillations were undetectable or only 
marginally detectable at times when the source exhibited 
clear episodes of dipping behavior, making it likely that 
the mechanism behind the dips causes the oscillations to 
become undetectable (Wijnands, Strohmayer & Franco 
2000). The dips can be understood as an obscuration 
of the inner system by matter periodically coming into 
the line of sight. This same matter then also attenuates 
the oscillations making them undetectable. Although no 
X-ray dips have been observed for both 4U 1728-34 and 
KS 1731-260, it is possible that the absence of NCOs at 
low mass accretion rates in these sources might also be 
due to changes in the environment. It is known that when 
M decreases in all atoll sources, the X-ray spectrum be- 
comes significantly harder and follows a power law spec- 
trum (see, e. g.. Barret et al. 2000, for a recent study of 
the hard X-ray spectrum of several atoll sources). It has 
been postulated that this power law spectrum is due to a 
comptonizing medium (e. g., a corona) around the inner 
system. If such a medium is indeed present at low M 
but not at high M, it might be responsible for smearing 
out the oscillations making them unobservable at low M . 
Similar mechanisms have also been introduced to explain 
the absence of the expected coherent pulsations in the 
persistent X-ray emission of neutron star LMXBs. 

However, this interpretation is not without problems. 
It has been shown that the only LMXB for which coher- 
ent millisecond oscillations have so far been found during 
the persistent emission (SAX J1808. 4-3658; Wijnands & 
van der Khs 1998) was in the island state (as judged both 
from the aperiodic rapid variability and the power law 
spectrum; Wijnands & van der Klis 1999; Gilfanov et al. 
1998) during the observations which led to its discovery as 
a millisecond X-ray pulsar. So, if a comptonizing medium 
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is present in 4U 1728-34 and KS 1731-260 which blurs the 
NCOs during the type-I X-ray bursts, why are the pulsa- 
tions in SAX J1808. 4-3658 not attenuated? This suggests 
that the mass accretion rate has a more direct influence 
on the production mechanism of the NCOs. 

Such a direct influence of M on the NCOs is also sug- 
gested by the fact that for seven bursts in 4U 1728-34 
the oscillations were visible during both the rise and the 
decay phases of the burst, while seven bursts show os- 
cillations only in the decay phase and two bursts show 
oscillations during the rise but not during the tail of the 
burst. Wc find that in 4U 1728-34 bursts with NCOs 
only during the tail occur at intermediate mass accretion 
rates, while bursts with NCOs in both the rise and the de- 
cay occur at high mass accretion rates. Our data further 
suggest that bursts with NCOs only in the rise occur at 
the highest mass accretion rates. Clearly, the exact mo- 
ment when the NCOs will be present during the bursts 
is determined by the mass accretion rate, which strongly 
suggests that the NCOs are indeed heavily influenced by 
the mass accretion rate directly and not by an M-induced 
environmental change. 

The behavior of the bursts, in particular the burst os- 
cillations, in 4U 1728-34 is very rich and many questions 
remain for which we have no adequate answers. Although 
the production mechanism behind the NCOs is not well 
understood, their coherence and their frequency stability 
over several years strongly suggest that indeed the oscil- 
lations are due to the spin frequency of the neutron star. 
But it is still difficult to imagine how the NCOs in the tail 
of bursts are formed. The entire surface of the neutron 
star is expected to be involved in the thermonuclear flash 
after the first few seconds. In such an environment a hot 
spot would probably not survive. A possible solution to 
this quandary is to invoke a different mechanism to pro- 
duce the NCOs in the decay phase of bursts. However, 
burst 17 in our sample (sec Fig. 2) shows what appears 
to be a continuous evolution of the NCO frequency. This 
strongly suggest that only one mechanism is behind the 
NCOs.' 

When we observe the oscillations in the tail of the burst, 
we usually also see an increase in the frequency by a few 
Hertz as the burst progresses (see also Strohmayer et al. 
1996; Strohmayer, Zhang & Swank 1998a; Strohmayer & 
Markwardt 1999). The observed evolution of the oscilla- 
tion frequency can be convincingly attributed to the ex- 
pansion of the neutron star atmosphere and conservation 
of angular momentum (see § 1). However, for at least one 
of our bursts (burst 17, see Fig. 2) we observe an increase 
in frequency of about 3.5 Hz. This would mean that the 
burning layer must have expanded by almost 50 meters. 
It is not clear if such a large expansion can be accounted 
for by the expanding layer interpretation, although the 
expansion can be at least as large as 40 meters (Cumming 
& Bildsten 2000). Such a large increase in frequency is 
not unique to 4U 1728-34. An increase of 5 Hz was re- 
cently found in the new NCOs source X 1658-298 (Wij- 
nands, Strohmayer & Franco 2000), which is even larger 
than what we have found for 4U 1728-34. However, be- 



cause of the higher frequency of the NCOs in X 1658-298 
(~567 Hz) compared to that in 4U 1728-34, the burning 
layer has to have expanded more in 4U 1728-34 than in 
X 1658-298 (50 meters versus 40-45 meters). 

When we investigated the strength of the oscillations 
in detail, we found that it is well correlated with inferred 
mass accretion rate: the oscillations increase in strength 
when the accretion rate of the system increases. This 
correlation is more evident within bursts observed only a 
few weeks apart, and becomes less clear for bursts sep- 
arated in time by several months to years. It is unclear 
why the correlation is strong on the short term (weeks) 
but fails to hold on longer timescales (months to years). 
Whatever the modulation mechanism, it will likely de- 
pend on the state of the fuel layers on the neutron star, 
which may change considerably with long term trends in 
the mass accretion rate. If wc think of the state of the 
neutron star surface as determining the moment for the 
onset of the modulation mechanism, then one can imagine 
that the overall changes in M over months and years are 
resurfacing the neutron star and hence the "critical" mo- 
ment will be reached at slightly different accretion rates 
on subsequent trips down the atoll track. In other words, 
if the time averaged accretion rate has been very low for 
a long interval of time, the structure of the neutron star 
surface might be different than when the time averaged 
accretion rate was considerably larger. 

The influences of the mass accretion rate on the burst 
oscillation properties strongly suggest that they are in- 
deed manifestations of physical properties of the atoll sys- 
tems. We would expect that this behavior of the NCOs 
will be also present in other sources. However, no detailed 
studies have been performed for the other sources similar 
to our study of 4U 1728-34 and that of KS 1731-260 by 
Muno et al. (2000). It is important to determine whether 
the same correlation will hold for all the sources which 
have NCOs. However, we can definitely exclude the pos- 
sibility that all X-ray burst sources will show NCOs when 
they are at high mass accretion rates. For at least two 
sources (Serpens X-1 and 4U 1705-44), X-ray bursts have 
been observed when they were on the banana branch, thus 
at relatively high mass accretion rates. No NCOs were 
discovered during those bursts (R. Wijnands 2000, private 
communication; Ford, van der Klis & Kaaret, 1998). A 
more intrinic reason (e. g., binary inclination, the strength 
and/or orientation of the magnetic field) must be found to 
explain the many burst sources that do not show NCOs. 
On the other hand, it is still possible that a considerable 
fraction of those sources have indeed only been observed 
in their island states and might still show NCOs when 
they can be observed to burst in their banana branches. 

To conclude it is clear that from the study of 4U 1728- 
34 and KS 1731-260 we can already glean some common 
behaviors and many differences. More atoll systems must 
be searched for NCOs in their bursts to begin differen- 
tiating between those trends that are universal accross 
atoll sources and those that are determined by intrinsic 
differences of the systems. 

Note added in manuscript. When we were about to 
submit our paper we became aware of the paper in prepa- 
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ration by van Straaten et al. (2000) presenting a similar 
analysis of 4U 1728-34. Although their study is limited 
to the AOl and A02 observations, they confirm our main 
conclusions: the NCOs are only present during bursts at 
the highest inferred mass accretion rates and the strength 
of the NCOs increases with accretion rate. They distin- 
guish between two types of PRE episodes, those marked 
by a small increase in radius after the main PRE episode, 
which they termed unusual, and the more standard con- 
tinuous radius decrease after the PRE episode. They 
find that standard PRE bursts only show NCOs after 
the peak of the burst and unusual PRE episodes occur 
only in bursts at the lowest M. We note that these two 
correlations also hold when including our A03 data set. 



I thank Craig Markwardt for extremely helpful and 
well documented software, Angela Olinto for unwavering 
support, Rudy Wijnands for several critical readings of 
the manuscript and Tod Strohmayer for unselfishly do- 
nating his data, time and knowledge and for teaching 
me all about RXTE and LMXBs. I thank Steve van 
Straaten and collaborators for their willingness to ex- 
change manuscripts before publication. This work made 
use of pubhc RXTE data obtained through the GSFC 
HEASARC archive, and was supported in part by NASA 
through grant NRA-99-01-ADP-195, by NSF through 
grant AST-0071235 by DOE through grant DE-FG0291. 



REFERENCES 



Arnaud, K.A., 1996, in ASP Conf. Ser. 101, Astronomical Data 

Analysis Software and Systems V, ed. G. Jacoby, & J. Barnes 

(San Francisco: ASP), 17 
Basinska, E. M., Lowin, W. H. G., Sztajno, M., Cominsky, L. R., 

& Marshall, F. J. 1984, ApJ, 281, 337 
Barret, D., Olive, J. F., Boirin, L., Done, C, Skinner, G. K., 

Grindlay, J. E. 2000, ApJ, 533, 329 
Bildsten, L. 1997, preprint (astro-ph/9709094) 
Bildsten, L. 2000, to appear in Proceedings of the 10th Annual 

October Astrophysics Conference, Cosmic Explosions, ed. S. S. 

Holt & W. W. Zhang, (New York:AIP), astro-ph/0001135 
Buccheri, R., et al. 1983, A&A, 128, 245 
Gumming, A. &: Bildsten, L. 2000, ApJ, submitted (astro- 

ph/0004347) 
Di Salvo, T., laria, R., Burderi, L., & Robba, N. R. 2000a, ApJ, in 

press (astro-ph/0005396) 
Di Salvo, T., Mendez, M., van der Klis, M., Ford, E., Robba, 

N. R. 2000b, ApJ, in press (astro-ph/0008361) 
Ford, E. G. & van der Khs, M. 1998, ApJ, 506, L39 
Ford, E. G., van der Klis, M., & Kaaret, P. 1998, ApJ, 498, L41 
Fujimoto, M., Y., Hanawa, T. & Miyaji, S. 1981, ApJ, 247, 267 
Gilfanov, M., Revnivtsev, M., Sunyaev, R. & Ghurazov, A&A, 338, 

L83 
Hanawa, T., & Fujimoto, M., Y. 1984, PASJ, 34, 495 
Hasinger, G., & van der Khs, M. 1989, A&A, 225, 79 
Lewin, H. G., van Paradijs, J. & Taam, R. 1993, Space Sci. Rev., 

62, 233 
Markwardt, G. B., Strohmayer, T. E., & Swank, J. H. 1999, ApJ, 

512, L125 
Mendez, M. 2000, in Nuc. Phys. B 80, XlXth Texas Symposium on 

Relativistic Astrophysics and Gosmology, ed. E. Aubourg, et al. 

(Netherlands: Elsevier) 15/16 GD-ROM (astro-ph/9903469) 
Mendez, M., & van der Klis, M. 1999, ApJ, 517, L51 
Mendez, M., van der Klis, M., Ford, E. G., Wijnands, R., & van 

Paradijs, J. 1999, ApJ, 511, L49 
Miller, J. M., Fox, D. W., DiMatteo, T., Wijnands, R., Belloni, 

T., Pooley D., Kouvehotou, G., & Lewin, W. H. G. 2000, ApJ, 

submitted (astro-ph/0003229) 



Miller, G. E. 1999, ApJ, 515, L77 

Muno, M. P., Fox, D. W., Morgan, E. H., & Bildsten, L. 2000, 

ApJ, submitted (astro-ph/0003229) 
Piraino, S., Santangelo, A., & Kaaret, P. 2000, A&A, in press 
Smith, D. A., Morgan, E. H., & Bradt, H. 1997, ApJ, 479, L137 
Strohmayer, T. E. 1999, ApJ, 523, L51 

Strohmayer, T. E. & Markwardt, G. B. 1999, ApJ, 516, L81 
Strohmayer, T. E., Zhang, W., Swank, J. H., Smale, A., Titarchuk, 

L., Day G. & Lee, U. 1996, ApJ, 469, L9 
Strohmayer, T. E., Jahoda, K., Giles, A. B., Lee, U., 1997a, ApJ, 

486, 355 
Strohmayer, T. E., Zhang, W., & Swank, J. H. 1997b, ApJ, 487, 

L77 
Strohmayer, T. E., Zhang, W., & Swank, J. H. 1998, ApJ, 503, 

L147 
Strohmayer, T. E., Zhang, W., Swank, J. H., White, N. E. & 

Lapidus, L 1998, ApJ, 498, L135 
Strohmayer, T. E. 1999, preprint (astro-ph/9911338) 
Swank, J. H., Becker, R. H., Boldt, E. A., Holt, S. S., Pravdo, S. 

H. & Serlemitsos, P. J. 1977, ApJ, 212, L73 
van der Klis, M. 1995 in X-ray Binaries, ed. W. H. G. Lewin, J. van 

Paradijs, & E. P. J. van den Heuvel, (New York: Gambridge 

University Press), 252 
van der Khs, M. 2000, ARAA, in press (astro-ph/0001167) 
van Paradijs, J. & Lewin, H. G. 1986, A&A, 157, LIO 
van Straaten, S., van der Klis, M., Kuulkers, E. & Mendez, M. 

2000, ApJ, submitted 
Wijnands, R., Strohmayer, T. E. & Franco, L. M. 2000, ApJL, 

submitted (astro-ph/0008526) 
Wijnands, R. & van der Klis, M. Nature, 394, 344 
Wijnands, R. & van der Klis, M. ApJ, 507, L63 
Zhang, W., Lapidus, I., Swank, L. H., White, N. E., Titarchuk, 

L. 1997, lAU Girc, 6541 
Zhang, W., Jahoda, K., Kelley, R. L., Strohmayer, T. E., Swank, 

J. H., Zhang, S. N. 1998, ApJ, 495, L9 



